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AN  OVERVIEW  OF  THE  COMMON  FLUID  MODELS  USED 
IN  FLUID-STRUCTURE  INTERACTIONS 


INTRODUCTION 


There  are  occasions  when  engineers  specialized  in  the  mechanics  of  solids  find  themselves  in  the 
position  of  having  to  model  the  response  of  a  solid  structure  which  is  interacting  with  an  adjacent  fluid 
whose  density  is  not  negligible  in  comparison  to  that  of  the  solid.  The  loading  on  the  solid  by  the  fluid  must 
then  be  properly  taken  into  account  in  order  to  accurately  predict  the  response  of  the  solid,  but  there  is  a  lack 
of  background  on  the  part  of  the  engineer  as  to  what  is  an  appropriate  fluid  model  to  couple  to  the  one  that 
he  is  already  using  for  the  solid.  The  overview  given  by  this  report  should  enable  such  an  engineer  to  obtain 
at  least  a  cursory  understanding  of  the  origin  and  limitations  of  the  fluid  models  which  are  in  common  use 
for  such  purposes.  With  this  as  a  base,  the  references  given  should  allow  one  to  pursue  both  greater  details 
and  actual  computer  code  implementations. 

This  report  derives  and  discusses  some  of  the  fluid  models  more  commonly  used  in  fluid-structure 
interaction  problems.  The  emphasis  here  is  on  linearized  equations  governing  disturbances  generated  in,  and 
propagated  through,  the  fluid  since  this  forms  such  a  large  part  of  the  fluid-structure  interaction  literature. 
The  fluid  may  either  surround  the  solid,  be  contained  by  the  solid  in  an  internal  cavity,  or  both.  In  many  of 
the  cases  in  which  the  fluid  surrounds  the  solid  the  fluid  is  idealized  as  extending  to  infinity.  This  is  a  source 
of  difficulty  in  numerically  modeling  the  physical  system  in  such  cases  since  one  cannot  directly  extend  the 
computational  grid  to  infinity.  Accurate  truncation  to  a  finite  volume  of  fluid  would  create  external  boundary 
condition  problems.  A  common  approach  to  circumventing  this  problem  in  the  linearized  case  is  to  refor¬ 
mulate  the  fluid  model  into  self-contained  equations  on  the  solid-fluid  boundary  which  directly  govern  the 
fluid  pressure  response  there.  Several  versions  of  this  approach  are  discussed  in  this  report.  They  are  all 
based  on  a  parent  fluid  model  which  is  derived  and  discussed  next 

FLUID  MECHANICS  FORMULATION 


The  behavior  of  the  fluid  in  fluid-structure  interaction  problems  is  governed  by  the  mass  conserva¬ 
tion  equation 


du 

Dp+ps ij 


^  =  0 


(1) 


the  momentum  conservation  equation 


(2) 


and  the  energy  conservation  equation 


where  the  usual  Einstein  convention  of  summation  over  repeated  subscripts  is  assumed,  and  where 
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p  =  fluid  mass  density, 


«,  =  fluid  velocity  vector, 

Fi  =  body  force  per  unit  mass, 
Gy  =  stress  tensor, 
e  =  internal  energy  per  unit  mass 
qt  -  heat  flux  vector 


u  =  uiUi 


.2 

D  m  1 
/  =  time,  and 


£  + 


x,.  =  spatial  position  vector. 

The  ‘*018161181  derivative”  Z?  can  be  interpreted  as  a  time  derivative  following  the  motion  of  the  fluid  because 

4-JpQdV  =  jpDQdV 
u  v  v 


holds  for  0  any  extensive  property  of  the  fluid  and  V  the  volume  of  any  fluid  particle  of  fixed  mass  moving 
with  the  flow.  Taking  the  vector  inner  product  erf  ut  with  (2)  and  subtracting  the  results  from  (3)  gives  the 
alternative  energy  equation 

du.  dq. 

pfle  =  ’’uWjSZj-  <4> 

The  above  governing  equations  must  be  supplemented  by  appropriate  equations -of-state,  by  appro¬ 
priate  constitutive  relations,  and  by  appropriate  boundary  conditions.  The  most  important  one  at  the  solid- 
fluid  interface  consists  of  the  no-slip,  no-penetration  boundary  condition 


«,  =  («*,) 


solid 


at  each  point  of  the  interface  at  each  instant  of  time.  The  stress  tensor  takes  the  form 


(5) 


Gy  =  -pby+dy 


(6) 


for  isotropic  fluids,  where 


p  *  (thermodynamic)  fluid  pressure  and 
dy  *  “viscous”  part  of  the  stress  tensor. 
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For  Newtonian  fluids  the  “viscous”  part  of  the  stress  tensor  is  linearly  related  to  the  spatial  derivatives  of 
the  fluid  velocity  as 


where 


Buk  f du:  dUj  2 


p.  =  (shear)  viscosity  and 

T)  =  expansion  (or  bulk)  viscosity. 


Substituting  (6)  into  (4)  and  then  using  (1)  leads  to  the 


version  of  the  energy  equation. 


The  nonequilibrium  setting  generally  present  in  a  flowing  fluid  requires  appropriate  interpretations 
[l,pp.  151-152]  for  thermodynamic  variables  such  as  entropy.  With  this  understanding,  the  thermodynamic 
relations 


Dp  =  c2Dp+  (|B)  £)s 

and 

TDs  =  De+pD(i) 

are  available,  for  which 

2  Bp 

c  =  s  local  speed  of  sound  squared, 

s  =  entropy  per  unit  mass,  and 


(8) 

(9) 


T  -  absolute  temperature. 


Equation  (9)  is  the  well  known  Gibbs  relation  for  a  simple  compressible  substance.  The  thermodynamic  re 
lation  [1,  p.  170] 

,Bpx  _  p c2pT 

where 


P  =  coefficient  of  thermal  expansion  and 
cp  s  specific  heat  (at  constant  pressure), 
combines  with  (8)  to  give 
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(10) 


Dp  =  c  Dp 


Combining  (1)  and  (10)  leads  to 


Combining  (2)  and  (6)  leads  to 


Combining  (7)  and  (9)  leads  to 


—2  pBT 

c  Dp*(’^j  =  -^Ds 


P  DUi 


dp  _  ddij 
<  +  =  pF,  +  ^y. 


3 xrv  1 


Bu:  Bq , 
pTDs  =  dl/s^-^. 


Equations  (11),  (12),  and  (13),  which  govern  the  fluid  behavior,  are  still  quite  general. 

The  first  assumption  usually  made  for  fluid-structure  interaction  problems  is  that  dissipative  and 
diffusive-type  processes  (of  molecular  transport  origin),  which  lead  to  internal  irreversibilities,  are  negligi¬ 
ble.  This  translates  into  the  specific  conditions 

9i  =  0 


dy*  0. 

which  represent  an  inviscid-fluid,  negligible-heat-conduction  postulate.  Implementing  these  conditions  in 
(13)  leads  to 

Ds  =  0  .  (14] 

Using  (14)  in  (11)  and  the  inviscid-fluid  condition  dy  —  0  in  (12)  gives 

Dp  +  pc2^  =  0  (15] 


pd“<+!  -  pF> 


respectively.  Equations  (15)  and  (16)  represent  the  governing  equations  for  the  fluid  behavior  under  the  giv¬ 
en  assumptions.  One  problem  which  originates  from  the  in  viscid  assumption  is  that  (5)  can  no  longer  be 
satisfied  in  its  entirety  in  conjunction  with  (15)  and  (16).  This  issue  will  be  addressed  more  fully  later. 

The  fluid  behavior  in  fluid-structure  interaction  problems  is  most  often  thought  of  in  terms  of  dis¬ 
turbances  in  the  fluid  which  propagate  through  it  The  fluid  models  can  be  categorized  according  to  whether 
the  disturbances  are  of  small  amplitude,  as  is  the  case  to  be  considered  in  this  report,  or  of  finite  amplitude 
[2, 3].  An  additional  consideration  is  whether  the  disturbances  are  produced  by,  or  interact  significantly 
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with,  a  background  flow  field  through  which  they  propagate  [4].  The  widely-used  fluid  model  emphasized 
in  this  report  assumes  that  the  total  flow  field  decomposes  into  a  steady,  mean  (background)  flow  field  of 
slowly  varying  velocity  and  small-amplitude  disturbances  superimposed  on  this  mean  flow.  This  decompo¬ 
sition  takes  the  form 


P  =  Po  +  Pl 

07) 

P  *  Po  +  Pl 

(18) 

ui  *  uiO  +  uil 

(19) 

2  2  2 
c*  =  ci+c\. 

(20) 

where  those  variables  labeled  with  a  0  subscript  are  associated  with  the  mean  flow  and  those  labeled  with  a 
1  subscript  are  associated  with  the  small-amplitude  disturbances  and  hence  they  have  small  magnitudes.  The 
steady  mean  flow  assumption  can  be  expressed  as 

dF,- 

5?  -  0  <21> 


a«.o 

Tt 


=  0 


(22) 


dpo 


~  0 , 


(23) 


whereas  the  slowly  varying  velocity  assumption  for  the  mean  flow  takes  the  form 


dui0 

57, 


=  0, 


(24) 


at  least  to  first  order.  The  mean  flow  is  independent  of  the  disturbances  and,  in  particular,  it  exists  even  in 
their  absence.  This  means  that  the  mean  flow  field  must  separately  satisfy  (13)  and  (16).  Substituting  (21) 
through  (24)  into  these  equations  leads  to 

=  0 
and 


3Po 

UJ0^j 


(25) 


These  two  equations  combine  to  give  the  relation 


PoMi0^i  =  0  • 


(26) 


Inserting  (17)  through  (20)  into  (13)  and  (16),  neglecting  second  order  terms,  and  utilizing  (22)  through  (26) 
leads  to  the  linearized  equations 
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(27) 


23l<n 

DoPi +Poeo5^ 


=  -Pt,u/iFi 


dpi 

po°o“n+5^ 


=  P/i  • 


where 


D» _ 

is  the  material  derivative  for  the  mean  flow. 

The  next  simplifying  assumption  consists  of  the 

~  «  1  (29] 

c0 

constraint  on  die  magnitude  of  the  body  force,  where 

F  -  a  representative  magnitude  of  the  body  force  per  unit  mass  and 
L  =  a  representative  length  scale  of  the  structure. 

The  constraint  (29)  is  certainly  true  for  the  most  common  case  for  which  the  body  force  is  gravity  [  1 ,  p.  1 6b] . 
The  pressure  gradient  term  in  (28)  is  of  the  order  [5,  p.  245] 

/  2^  \ 


pp.l-ofr1) 


in  magnitude  and  hence  the  estimate 


WpA 


■1  UJ 


follows  for  equation  (28).  (The  II A  |j  symbol  represents  the  norm  of  A  for  any  given  A  with  components 
Ar  In  this  case  it  is  the  Euclidean  norm  so  that  it  represents  the  magnitude  of  the  vector.)  The  estimate 

Po“;|P/  cfLF\ 

is  also  obtained,  for  equation  (27),  so  that  (27)  and  (28)  reduce  to 

7dun 

DoPi  +  Poco5Ty  =  0  (30 


dpt 

Pt/Vto+gj"  -  0 
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under  the  constraint  of  (29).  Equations  (30)  and  (3 1),  which  govern  the  fluid  behavior  under  all  of  the  above 
assumptions,  are  valid  for  an  inhomogeneous  fluid  for  which  pQ  and  c\  are  functions  of  position. 

Equation  (31)  can  be  integrated  easily  for  the  common  case  in  which  pQ  is  taken  to  be  a  constant 
and  the  gradient  of  a  velocity  potential  is  utilized  for  uiX .  This  assumption,  along  with  the  assumption 
of  constant  Cq,  will  be  taken  as  valid  for  the  remainder  of  the  report.  This  leads  to 


=/('> 

as  the  integration  of  (31)  for  an  arbitrary  function /(/)  of  time  only.  There  is  no  loss  in  generality  in  taking 
/  to  be  zero  since  any  nonzero  /  can  be  absorbed  into  the  definition  of  4>j  as 

without  affecting  uiX .  The  above  expression  hence  reduces  to 


Pi  =  -Po0O*i 


(32) 


with  the  velocity  potential  defined  by 


=  5Z  • 


(33) 


Equation  (32)  can  be  substituted  into  (30)  to  give 

2  2^2^i 

°»*'-c»3^;  =  0  (34) 

which,  with  appropriate  boundary  conditions,  can  be  solved  for  the  disturbance  solution  <{>, .  The  complete 
model  for  the  fluid  disturbance  is  hence  given,  at  this  point,  by  equations  (32)  through  (34)  and  appropriate 
boundary  conditions.  Operating  on  (34)  with  -p0OQ  and  using  (32)  leads  to 


^>-4^=°  (3s 

as  the  equation  for  px . 

Classical  Wave  Equation  of  Acoustics 

When  the  above  equations  are  transformed  into  a  reference  frame  which  is  at  rest  with  respect  to 
the  fluid  they  then  describe  the  fluid  as  an  acoustical  medium  whose  behavior  is  governed  by  the  classical 
wave  equation.  This  is  the  simplest,  and  probably  the  most  common,  rigorous  fluid  model  in  use  for  fluid- 
structure  interaction  problems.  It,  and  approximations  to  it,  will  form  the  basis  of  the  remainder  of  the  report. 
The  simplest  way  to  obtain  the  version  of  the  above  fluid  model  in  a  frame  at  rest  with  respect  to  the  fluid 
is  to  take  ui0  =  0  in  all  of  the  above  expressions  so  that  DQ  reduces  to  d/dr.  To  ease  the  notational  burden, 
the  subscript  1  on  pv  $j,  and  uiX  will  be  dropped.  In  addition,  since  pj  and  cx  are  never  referred  to,  the 
subscript  0  on  pQ  and  c0  will  be  dropped  so  that  p  will  represent  the  constant  p0  and  c  will  represent  the 
constant  c0.  The  p  of  (17)  will  be  denoted  by  Ptotal  so  that  (17)  becomes 
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P  total  ~  Po+P- 


(36) 


Equations  (23)  and  (25)  are  repeated  here  as 


dPo 

Tt 


=  0. 


(37) 


and 


3P°  - 

as;  ~  pF‘ 


(38) 


under  this  new  notation.  (Equation  (21)  is  not  needed  since  it  is  implied  by  (37),  (38),  and  die  fact  that  p  is 
a  positive  constant)  Equations  (32),  (33),  and  (34),  which  govern  die  fluid  disturbance  propagation,  become 


p  =  ‘p5 1 


(39) 


“«  = 


5V 


(40) 


and 


924> 

a?~c  Zxfc 


=  0. 


As  in  the  case  for  4>,  the  classical  wave  equation 


if 


c\ ^  "  0  ’ 


(41) 


(42) 


results,  in  this  case  from  (35),  for  the  acoustic  pressure  p .  The  disturbance  part  of  the  total  mass  density  also 
obeys  the  same  wave  equation,  but  it  is  not  usually  explicitly  used  in  modeling  the  fluid. 


Other  than  boundary  conditions,  which  will  be  discussed  in  the  next  subsection,  equations  (36) 
through  (41)  give  a  complete  model,  henceforth  referred  to  as  the  potential-based  model,  of  the  fluid  behav¬ 
ior.  Equation  (42)  is  an  equivalent  substitution  for  (41)  in  the  above  complete  set,  for  t  >  tQ,  if  the  initial 
condition 


la/2  C  dxjdxj  Jt 


=  0 


*0 


(43) 


is  imposed  on  This  is  most  often  satisfied  by  the  condition  that  0  is  a  constant,  usually  zero,  at  some  time 
/  =  /0  and  earlier.  An  alternative  fluid  model,  henceforth  called  the  pressure-based  model,  is  hence  given 
by  equations  (36)  through  (40),  (42),  and  (43).  The  solution  to  the  governing  equations  ultimately  provides 
a  history  of  the  distribution  of  Ptotai  along  the  solid-fluid  interface,  which  is  required  for  predicting  the  re¬ 
sponse  of  the  solid. 
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Boundary  Conditions 

Either  of  the  above  alternative,  equivalent  models  is  incomplete  without  appropriate  boundary  con¬ 
ditions.  The  boundary  condition  for  p0  is  simply 

P0  (xio)  ~  Poo  •  (44) 

where 


xi0  =  a  given,  fixed  point  in  the  fluid  or  on  its  boundary  and 
p0 o  =  a  prescribed  constant  value. 

Both  of  the  alternative  models  are  based  on  solving  the  classical  wave  equation 

dV  2^V 


a  t2  c 


=  0 , 


(45) 


where  V  is  either  <j>  or  p,  depending  on  the  model.  The  appropriate  boundary  information  for  this  equation 
is  to  specify  either  y  (Dirichlet),  dy/dn  (Neumann),  or  a  linear  combination  (mixed)  of  them  on  the 
boundary,  where 

=  and 

ni  =  the  unit  normal  vector  on  the  boundary,  specified  as  either  inward  or  (usually)  outward. 

In  addition  to  this,  one  needs  to  specify  an  initial  condition  for  y  (jr(-,  t)  as 

y  (*,-,  0)  =  y0(x,)  ,  (46) 

where 

y0  (Xj)  is  the  prescribed  intial  distribution  of  y . 

In  each  model  the  closest  approximation  that  one  can  make  to  (5)  at  the  solid-fluid  interface  requires 
the  Neumann  boundary  condition.  Taking  the  ordinary  vector  inner  product  of  (5)  with  n,  gives  the  result 


u-  ■  <“->„«• 


(47) 


where 


un  -  tt jUj  is  defined  as  the  normal  component  of  the  velocity  vector. 

Iklring  d/d/  of  (5)  and  then  taking  the  vector  inner  product  of  the  results  with  ni  leads  to 

an  ~  (°n)  solid' 


(48) 
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where 


Bu 


an  =  n 


is  defined  as  die  norma]  component  of  the  acceleration  vector. 


Both  (47)  and  (48)  are  only  approximations  to  (5)  since  they  each  contain  only  part  of  the  information  con¬ 
tent  of  (S).  Equations  (47)  and  (48)  are  not  equivalent  as  can  be  seen  from  the  identity 

Bur  Bn, 

57 


WH  '"V  . 


(49) 


The  vector  Bn  fit  is  easily  seen  to  be  tangential  because 

Bn, 


n  ni .  b  (njnj\  _  _  o 

nJTt  "57(^")  ~  5r(2)  "°- 


If  Bn  fit  is  negligible  or  zero  then  taking  B/Bt  of  (47)  gives  (48),  and  integrating  (48)  over  time  from  tQ 
to  t  gives  (47)  for  t  >  t0  if  (47)  is  true  at  t  =  tQ.  This  is  at  least  approximately  the  case  when  the  displace¬ 
ments  of  the  solid  during  its  response  are  small.  The  relation  (47)  is  used  for  the  potential-based  fluid  model 
as  can  be  seen  by  taking  the  vector  inner  product  of  (40)  with  nt  and  substituting  (47)  into  the  results  to  get 

si  =  <*» 

as  the  appropriate  fluid-solid  interface  boundary  condition.  Taking  the  gradient  of  (39),  taking  the  vector 
inner  product  of  the  results  with  nt,  and  substituting  (48)  into  the  results  leads  to 

%  *  sou  (51> 

as  the  appropriate  fluid-solid  interface  boundary  condition  for  the  pressure-based  fluid  model.  The  two  mod¬ 
els  hence  lead  to  different  results  when  the  unit  normal  at  the  interface  is  not  constant  in  time. 


In  many  of  the  problems  in  fluid-structure  interactions  the  fluid  surrounds  the  solid  and  goes  off  to 
infinity  in  all  directions.  In  these  cases  one  needs  to  consider  an  additional  boundary  condition  which  applies 
as  one  approaches  infinity.  A  radiation  boundary  condition,  which  will  be  discussed  later  where  its  meaning 
will  be  easier  to  understand,  is  usually  appropriate. 

REFORMULATION  INTO  INTERFACE  INTEGRAL  EQUATION 


The  source  of  the  disturbances  propagating  through  the  fluid  may  be  the  solid  itself,  or  the  source 
may  be  located  within  the  fluid,  or  both.  The  modeling  of  incident  pressure  waves  which  are  emanating  from 
a  source  within  the  fluid  and  impinging  on  the  solid  can  be  handled  by  including  a  source  term  O  (of  bound¬ 
ed  extent)  in  the  governing  wave  equation  (45)  for  the  disturbances,  so  that 


aV  -2  a  V  _  „ 

a x£xj~c 


(52) 


governs  the  disturbances  for  y  either  the  velocity  potential  or  the  acoustic  pressure,  where  0  =  0  through¬ 
out  most  of  the  fluid.  Let  y  denote  the  Fourier  transform,  with  respect  to  time,  of  y ,  so  that 


y(<o)  =  (27t)-v2j!^.y(0«/>(»0)0d/ 


to 


and  similarly  for  any  other  dependent  variables.  As  in  [6],  equation  (52)  will  be  reformulated  into  a  surface 
integral  equation  through  the  intermediate  use  of  Fourier  transforms.  Diking  the  transform  of  (52)  and  using 

gj*  =  -t©V((D)  (53) 


leads  to 


3^  ,2- 


where 


k  ~  to/ c  defines  k. 


An  equation  governing  the  solution  to  (54)  on  the  surface  bounding  the  fluid  can  be  obtained 
through  the  use  of  a  Green's  function  G  defined  by 


where 


(5^+*2)c(x'lX'<)  =  «<*.-*'(>  • 


8  *  the  Dirac  delta  function,  (actually  a  generalized  function  or  distribution). 


For  the  case  where  (55)  is  solved  for  an  infinite  space  one  gets  the  solution 


where  rt  and  r  are  defined  by 


/■*,_  _  exp(ikr) 

G(x>x‘)  s — 4 ST" 


ri  =  xi-x'i 

r  =  ('V*)1'2- 


The  two  identities 


1 IIC  IWU  1UCUUUC5 

and 

which  follows  from  (54)  and  (55),  can  be  combined  to  give 

yS  =  C5  +  ^(v|)  .  (59, 

Let  SR  denote  the  surface  of  a  sphere  of  radius  R  whose  fixed  center  lies  within  the  solid.  The  radius  R  is 
taken  to  be  large  (ultimately  infinity)  so  that  the  surface  SR  is  completely  within  the  fluid  and  it  completely 


u 


contains  all  of  the  source  C.  (The  source  <5  is  zero  outside  of  SR.)  Let  S  denote  the  solid-fluid  interface 
surface  which  encloses  the  solid.  Let  ST  denote  the  surface  defined  by 


ST  =  S\jSr 

and  denote  the  fluid  volume  which  ST  encloses  by  V.  Integrating  (39)  over  V  leads  to 

SdY+j  (y^L-G^dS' ,  (60) 

where  the  unit  normal  in  the  surface  integral  is  outward  with  repect  to  the  fluid,  upon  using  Gauss’  theorem, 
which  states  that 


fySc/V  =  f 
v  v 


v  i  s 

for  reasonably-behaved,  but  otherwise  arbitrary^,-  and  for  any  (reasonable)  given  volume  V  which  is  en¬ 
closed  by  surface  S  with  outward  unit  normal  n(-.  Let  VSource  denote  the  volume  of  the  source  a,  so  that 
O  is  only  nonzero  within  Vsource-  Splitting  ST  into  its  component  surfaces  SR  and  S,  equation  (60)  be¬ 
comes 


JvSdV  =  j  GodV+J(v||-C^)dS'+f()ji||-G^)<«'.  (61) 

^  V»mrct  &  $K 

The  radiation  boundary  condition,  with  G  from  (56),  leads  to  the  vanishing  of  the  SR  integral  in 
(61)  as  R  approaches  infinity.  One  formulation  of  the  radiation  boundary  condition  is  to  take  the  \jr  on  SR 
to  be  a  purely  outgoing  wave  of  the  form 


V  =  Vo*0’** 


exp(ikR) 

R 


(62) 


where  \yQ  satisfies 


^(sine^0) 


=  0 


so  that  V  satisfies  the  homogeneous  (O  equal  zero)  version  of  (54), 


Ri  *  x'i~xiO' 

R  =  {RjRj)m . 
x'i  is  on  SR 

xiQ  is  the  center  of  the  sphere  SR  ,  and 

8  and  ^  are  the  spherical  polar  angles  specifying  the  direction  of  . 

(The  equation  for  \|/Q  can  be  transformed  into  Laplace's  equation  in  p  and  $  by  the  substitution  =  log[(l- 
cos8)/sin0].)  For  large  R  one  finds  that 


and 


to  first  order  in  t/R ,  where 


ei  ®  xi  ~  xi0  • 

e  =  (eft) 1/2  , 

rt  =  £.  -  Rj  for  the  rf  of  (57), 

jc4-  is  on  S  , 


-(w>”-jt[i  +  (Ljfi)]  , 


s*1  • 

These  approximations  to  r  and  r-1  are  substituted  into  (56)  and  the  resulting  G,  along  with  (62),  are  then 
substituted  into  the  SR  integral  term  of  (61).  It  is  easily  seen  that  the  lowest  order  (r  **  R)  approximation 
gives  identically  zero  for  this  SR  integral  term  and  hence  the  first  order  terms  must  be  retained.  The  result 
is  that  the  SR  integral  in  (61)  is  of  order  MR  so  that  the  integral  vanishes  as  R  goes  to  infinity. 

The  a  term  in  (61)  represents  the  waves  emanating  from  the  (transformed)  source  o  in  an  infinite- 
space  fluid,  that  is,  in  the  absence  of  the  solid.  As  such,  it  represents  the  incident  part  of  the  wave  which 
impinges  on  the  solid  from  sources  internal  to  the  fluid.  It  is  hence  relabeled  as 


v,„=  J  GadY 


where  the  subscript  inc  denotes  the  incident  part  of  the  wave.  Taking  the  inverse  transform  of  (63)  and  us¬ 
ing  the  relation 

f(t-a)  =  exp(i(oa)f((D)  (64) 

for  arbitrary  (reasonably  behaved)  /,  with  /  =  O  and  a  =  r/c  since  k  =  to/c,  leads  to  the  well  known 
infinite  space  result 

1  .  o  (*'.•,  t-r/c) 

v*  (*,.<)  -  -j;  J  - - - dV  (65) 

^mrei 

upon  using  (56),  where  r  is  given  by  (58).  The  yj||C  of  (65)  satisfies  (52)  in  an  infinite  fluid  (no  solid). 

When  the  radiation  boundary  condition  and  (63)  are  inserted  into  (61)  and  R  is  then  taken  to  infinity 
one  gets  the  result 
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Substituting  (56)  into  this  result  gives 

JV&ZV’  =  yf(|c+  i-J  (r'V/yCitr-lJijf+r'1^}  exp(ikr)dS' 

upon  using 

dr  n'fj 
M  =  r  ’ 


(66) 


where  ri  and  r  are  defined  by  (57)  and  (58)  respectively.  Switching  to  an  inward  unit  normal  vector  with 
respect  to  the  fluid  gives 

f V&/V  =  \flHC  +^J  {  r^n'jtj  ( 1  -  ikr)  y  -  r'1^}  exp  ( ikr )  dS' ,  (67) 

r  S 

where  the  n\  vector  is  now  outward  with  respect  to  the  solid.  The  delta  function,  at  least  the  specific  one 
defined  by  (55)  through  (58),  has  the  property  that 


/(*,.)  for  xi  inside  V 
=  •  0  fOT  x> ou,sid' v 

4i/(Xi)  for  X‘  “  S 


(68) 


for  continuous  functions  /,  where  S  is  the  surface  enclosing  the  given  volume  V  and 

£2  =  (the  solid  angle  subtended  by  the  surface  5  at  JCf)  =  2 it  for  smooth  surfaces. 

(The  first  two  properties  are  universal  for  all  delta  functions.)  Using  V  for /in  (68)  and  substituting  the 
results  into  (67)  gives 

V(x,,*)  =  2yinc  +  {  r-’n'fj  ( 1  -  ikr)  V  "  exp  (ikr)  dS'  (69) 

as  the  governing  integral  equation  for  y  on  S  assuming  that  S  is  smooth, 

V(*,,0  =  VIIIC  +  ^  J  {r~Vyr.  ( 1 "  ikr'*  V  ~  exp  (ikr)  dS'  (70) 

as  the  equation  for  finding  y  within  the  fluid  after  having  first  found  it  on  S  from  solving  (69),  and 

0  *  < 1  ~ ikr)  V  -  exP  (**r)  *8'  (?1) 

s 

as  the  equation  for  xi  within  the  solid.  Equation  (71)  serves  as  an  alternative  to  (69)  for  determining  y  on 
S.  The  transform  of  (50)  with  y  =  $  for  the  potential-based  model,  or  the  transform  of  (51)  with  y  «  p  for 
the  pressure-based  model,  are  substituted  for  dy/dn'  in  (69)  to  obtain  a  governing  integral  equation  for  the 
fluid  response  on  5  which  includes  the  boudary  conditions. 
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Solving  the  Equations  in  Frequency  Space 

In  some  of  the  cases  for  which  either  the  solid  constitutive  behavior  is  linear  or  the  solid  behavior 
is  prescribed  it  is  useful  to  solve  for  the  fluid  response  in  the  frequency  space  setting  represented  by  (69) 
and  the  transform  of  either  (50)  or  (51).  It  is  well  known  [7],  however,  that  (69)  with  prescribed  1/dn’ 
and  zero  y.  does  not  yield  a  unique  solution  at  certain  discrete  values  of  k  (that  is,  characteristic  frequen¬ 
cies,  since  k  =  co/c).  Reference  [7]  offers  the  remedy  of  using  a  discretized  version  of  (71)  at  selected 
locations  to  supplement  the  main  model,  a  discretized  version  of  (69),  so  as  to  form  an  overdetermined  sys¬ 
tem  which  is  solved  by  least  squares.  This  approach, which  is  implemented  in  the  computer  code  CHIEF,  is 
able  to  provide  accurate  solutions  even  at  the  problem  frequencies.  Alternative  approaches  exist  [8  and  9] 
for  overcoming  the  same  problem,  such  as  that  embodied  in  the  computer  code  CONDOR. 

Time  Domain  (Retarded)  Formulation 


Equation  (69)  can  be  reformulated  into 

v<*,*>  =  2V(.c+^J{r-V/y(v+t|')-r-,||)«p(,*r)1iS’ 

s 

using  (53).  The  inverse  transform  of  this  equation  gives 

y  (xif  t)  =  2yi||C  +  ^  J  { r  3 n'fj  ( <V>  +  ~  (gf))  “  r  1  (72> 

upon  using  (64),  where  the  pointed  brackets  denote  evaluation  at  the  retarded  time  as 

(f)  =  f  (x'it  t  —  -  )  for  an  arbitrary  function  /  of  position  and  time, 

C 

r  is  given  by  (57)  and  (58),  and 


n\  is  the  outward  unit  normal  with  respect  to  the  solid. 
Equation  (72)  can  be  combined  with  (50)  to  produce 

a  /  _  1  f  r  - 3_«  _  /  /a\  i  ^  /^\\  n  jo 


♦  <*,'>  =  2«w  +  2*1  <r"V/>  «*>  +  J  <S  »  >  ^  -  2il  ( ^  <73> 

5  S 

for  y  =  4>  in  the  potential-based  model,  or  with  (51)  to  produce 

P  (xf  0  =  +  +  <74) 

s  s 

for  y  =  p  in  the  pressure-based  model.  Methods  which  are  based  on  the  solution  of  either  (73)  or  (74)  are 
referred  to  as  retarded  potential  methods.  Computer  code  descriptions  based  on  discretized  versions  of  either 
(73)  or  (74)  can  be  found  in  references  [10-12]. 


APPROXIMATIONS  TO  RETARDED  POTENTIAL  METHOD 


There  are  three  common  approximations  to  the  retarded  potential  method  which  are  referred  to  as 
the  early-time  approximation,  the  late-time  approximation,  and  the  DAA  (Doubly  Asymptotic  Approxima¬ 
tion). 
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Early-Time  Approximation 

If  S  is  planar  and  pinc  is  taken  to  be  zero  then  (74)  reduces  to 

P(x„0  <tS'  (75) 

since  n'  -Tj  -  0  everywhere,  where  n\  is  now  a  constant  unit  vector  along  S  and  S  =  (plane  of  z-  con¬ 
stant).  If,  starting  from  zero,  t  represents  the  elapsed  time  and  if  tire  solid  surface  normal  acceleration  is  zero 
for  negative  times  then  the  response  at  any  given  point  on  S  can  be  influenced  only  by  points  of  5  within  a 
radius  of  ct.  Equation  (75)  hence  reduces  further  to 

p  (*, y,  t)  =  ^/tyjo  JJ*fly(x  +  rcos0,y  +  rsin0,/-^)dedr  (76) 

upon  using  dS'  =  rdrdQ  with  planar  polar  coordinates,  where  it  is  understood  that  a,  represents  the  solid 
surface  acceleration.  The  early-time  approximation  assumes  that 

ct«L 

where L  is  a  characteristic  length  scale  for  5.  The  fact  that  0  <.  r  £  cf  bounds  r  leads  to  r  «  L,  so  that  a 
Taylor  series  expansion  gives 

a.-  (x  +  rcos0,  y  +  rsinO,  t  - -)  «  a,  ( x ,  y,  t  - -)  +  O  (r) 

c  1  c 

for  the  ai  in  (76).  (Note  that  r/c  is  not  small  compared  to  /  in  this  case  and  hence  no  expansion  is  warranted 
for  the  time  dependence.)  For  small  enough  t,  the  radius  Ct  becomes  small  enough  that  the  initial  assump¬ 
tion  of  a  planar  surface  becomes  true  locally  even  for  an  5  with  global  curvature,  at  least  as  long  as  the  sur¬ 
face  is  smooth.  Equation  (76)  now  reduces  to 

P  (*. y,  t)  =  y,t-r-)dr 

since  there  is  no  longer  any  0  dependence  in  the  integrand.  The  change  of  variable  t'  =  t  —  r/c  transforms 
the  above  integral  into 


P  (■*.  y,  0  =  Pcnjjoaj  (*’  y,  O  dt'  >  (77) 

which  is  one  form  of  the  early-time  approximation.  Differentiating  (77)  with  respect  to  time  gives  another 
useful  form 

I  =  P'KU,  <78> 

on  5  for  the  early-time  approximation.  Higher  order  results  for  small  ct  can  be  found  in  reference  [13]. 
Late-Time  Approximation 

The  late-time  approximation  is  based  on  the  assumption 

ct»L , 
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which  leads  to  t»  r/c  because  r  is  ultimately  bound  by  r  £  L  in  magnitude  if  L  is  the  length  scale  asso¬ 
ciated  with  the  largest  global  extent  of  S.  (The  above  approximation  can  also  be  thought  of  in  terms  of  the 
wavelength  X  so  that  c/(0  =  X  »  L.)  A  Thylor  expansion  for  small  r/c  for  generic  retarded  functions  f 

<f>  i)  -  +  |r<xV'>  +0<£)3  09) 

is  hence  relevant  in  this  case.  Using  (79)  in  (72)  with  zero  for  yinc  gives 

_  If  r  3  _  t  _  ...  -l^V  i  mi  snn\ 


upon  retaining  only  lowest  order  terms.  Equation  (80)  can  be  derived  from  Laplace’s  equation 

dV 

zfc  - 0  <81> 

in  direct  analogy  to  the  derivation  of  (72)  from  the  wave  equation  (52)  with  no  internal  sources  (so  that  yinc 
and  0  are  zero).  Equation  (80)  can  hence  be  thought  of  as  representing  the  governing  equation  on  5  for 
solutions  to  (81).  If  y  is  taken  to  be  the  velocity  potential  $  then  (80)  describes  a  potential  flow  and  the 
flow  is  hence  incompressible  near  the  solid  body  [5,  p.  280].  Equation  (81),  equation  (50)  on  S,  and  the  van¬ 
ishing  of  the  gradient  of  $  at  infinity  together  form  a  boundary  value  problem  for  $  which  can  be  reformu¬ 
lated  in  terms  of  a  distributed  surface  source  on  5  so  as  to  ultimately  produce  an  added  mass  model  [14]  of 
the  fluid  response.  An  alternative  route  to  an  added  mass  model  is  initiated  by  taking  y  =  p  in  (80)  to  get 


p  = 


where 


1  f  —3 

r  (/)  *  {r  n'jTjf]  dS'  for  T  operating  on  a  generic  function  /, 
P  f  “"1 

O  (/)  =  — J  (r  f)  dS'  for  <X>  operating  on  a  generic  function  f,  and 


r/  is  a  shorthand  notation  for  T  if)  and  similarly  for  G>. 

The  use  of  an  operator  notation,  where  the  integrations  of  (80)  are  symbolically  represented  as  the  operators 
r  and  O  in  (82),  is  common  in  functional  analysis  and  allows  a  great  economy  and  conciseness  in  the  de¬ 
velopment  Inserting  (51)  into  (82)  leads  to 


p  =  r 


which  can  be  symbolically  solved  for  p  as 


P  =  MO 


*'  solid 


where 


/  =  the  identity  operator  and 
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A  =  (/-D^Ofor  (/-T)”1  as  the  inverse  of  the  operator  (I-T). 
The  inverse  T~l  of  an  operator  T,  when  it  exists,  has  the  property 


rxTf  =  rrlf  =  / 

for  a  generic  function  /,  the  Laplace  transform  and  inverse  transform  being  examples  of  such  a  T  and  T~l 
respectively.  An  alternative,  less  abstract,  interpretation  consists  of  taking  T  and  O  as  discretized  versions 
of  die  cot  responding  integral  operators  of  (80)  so  that  T  and  <S>  are  then  matrices  and  p,  dp /Bn,  and  the 
solid  normal  surface  acceleration  are  then  column  vectors.  Under  such  an  interpretation,  which  would  also 
be  the  practical  approach  in  an  actual  implementation,  equation  (82)  becomes  an  approximation  to  (80),  the 
/  becomes  the  identity  matrix, 

(/ -  O'1  becomes  the  inverse  of  the  matrix  (I-T) ,  and 

A  is  the  matrix  multiplication  of  (/  -  T) -1  and  <D. 

The  presentation  in  [14]  is,  in  fact,  given  in  terms  of  matrices  originating  from  a  finite  element  discretiza¬ 
tion.  The  force  on  the  solid  per  unit  area  due  to  the  fluid  is  -pni  and  hence,  by  Newton’s  second  law  and 
(83),  the  A  can  be  interpreted  as  an  “added  mass  per  unit  area”  operator  (or  matrix)  which  is  due  to  the  pres¬ 
ence  of  the  fluid.  In  other  words,  for  an  element  at  the  surface  S  in  a  finite  element  subdivision  of  the  solid, 
adding  the  fluid  force  acting  on  this  element  to  the  other  such  forces  is  equivalent  to  adding  a  mass  operator 
to  the  mass  of  the  solid  for  that  element  in  die  equations  of  motion. 

Doubly  Asymptotic  Approximation 

The  Doubly  Asymptotic  Approximation  (DAA)  originated  with  Geers  [15]  and  represents  a  smooth 
joining  of  the  early-time  and  late-time  (in  the  form  of  an  added  mass  model)  approximations.  The  simplest 
such  joining  consists  of 

A^  +  pcp  =  pcA(aH)sol.d,  (84) 

so  that  (84)  reduces  to  (78)  when  the  dp/dt  term  dominates  and  (84)  reduces  to  (83)  when  the  p  term  dom¬ 
inates.  (It  is  assumed  that  the  inverse  of  A  exists.)  Equation  (84)  is  denoted  by  DAAj,  the  next  higher  order 
version  being  denoted  by  DAA2  [16],  An  advantage  of  the  DAA  over  the  retarded  potential  method  is  that 
the  DAA  equation  is  local  in  time.  This  leads  to  substantial  computational  savings  for  the  DAA  when  the 
methods  are  implemented  numerically.  The  tradeoff  is  that  the  DAA  is  less  accurate  than  the  retarded  po¬ 
tential  method  in  a  variety  of  circumstances  [17]. 
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